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SUMMARY 


OBJECTIVE 

Determine  the  probable  cause  for  a  discrepancy  in  the  Consultative  Committee 
International  Radio  (CCIR)  Report  322-3  radio  noise  model  and  recommend  a  course 
of  action  to  overcome  the  discrepancy. 

RESULTS 

The  basis  for  this  discrepancy  was  found  to  be  in  the  procedure  used  to  prepare 
the  measured  noise  data  for  the  determination  of  a  global  numerical  representation  of 
the  1-MHz  data.  The  procedure  followed  in  the  development  of  the  model  was  to 
determine  correction  factors  to  the  old  CCIR  model  for  each  measurement  site,  to 
interpolate  these  corrections  to  a  100-latitude  by  84-longitude  grid  for  each  time 
block/season,  to  add  the  correction  factors  at  each  grid  point  to  corresponding  values 
for  the  old  CCIR  model,  and  finally  to  numerically  map  the  resulting  data  for  each 
time  block  and  season.  Nineteen  locations  were  used  in  the  final  model.  Four  sites 
used  in  the  original  CCIR  model  were  not  used.  These  include  Bill,  Wyoming;  Byrd 
Station,  Antarctica;  Ibadan,  Nigeria;  and  Thule,  Greenland.  As  no  correction  factors 
were  obtained  for  these  locations  nor  a  correction  factor  of  zero  used,  the  interpolation 
algorithm  used  to  obtain  the  100-latitude  by  84-longitude  grid  of  correction  factors 
supplied  other  values.  For  Bill,  Wyoming,  the  result  is  not  too  serious;  but  for  the 
other  three  sites,  the  error  is  at  some  seasons  and  time  of  day  serious.  For  Thule, 
Greenland,  the  maximum  and  minimum  errors  in  the  correction  contours  were  10.1 
and  -10.8  dB,  respectively.  For  Ibadan,  Nigeria,  the  maximum  and  minimum  errors 
were  12.5  and  -1.5  dB,  respectively.  For  Byrd  Station,  Antarctica,  the  maximum  and 
minimum  errors  were  12.0  and  3.0  dB,  respectively.  Examination  of  the  geographical 
extent  of  these  errors  reveals  that  the  error  is  not  confined  to  the  measurement 
location  but  in  fact  is  very  large.  It  was  found  that  the  error  as  a  function  of  frequency 
was  diumally  dependent.  An  error  of  10  dB  at  1  MHz  was  more  serious  at  another 
frequency  during  local  daytime  than  at  night.  Finally,  the  absence  of  the  data  locations 
affected  the  accuracy  of  the  interpolation  itself. 

RECOMMENDATIONS 

1.  Use  the  CCIR  Report  322-3  atmospheric  noise  model  with  caution,  especially  in 
the  northern  and  southern  high  latitudes,  the  Arabian  Peninsula,  northern  Africa,  and 
the  mid-Atlantic  Ocean  area.  In  these  areas,  consider  using  the  original  CCIR  Report 
322  model. 

2.  A  three-step  process  should  be  followed  to  develop  a  new  1-MHz  atmospheric 
noise  model. 

3.  First,  obtain  correction  factors  for  additional  locations  to  increase  the  accuracy 
of  the  interpolation. 


i 
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4.  Second,  test  the  method  of  interpolation  against  a  benchmark. 

5.  Third,  use  the  2Lacharisen  and  Jones  numerical  mapping  technique  applied  in 
local  time  to  develop  the  final  model. 

6.  Consider  using  a  latitude  transformation  to  increase  the  accuracy  of  the 
numerical  mapping  technique. 

7.  Submit  a  corrected  model  to  the  CCIR. 
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1.0  INTRODUCTION 


Recently,  Bowen  and  Fraser-Smith  (1992)  made  a  comparison  of  measured  32-kHz 
radio  noise  amplitudes  with  the  current  CCER.  Report  322-3  noise  model  predictions 
(CCIR,  1988).  They  found  that  the  greatest  discrepancies  between  the  measured  and 
predicted  amplitudes  were  observed  at  the  two  northern  high-latitude  stations  (Sdndre 
Str^fjord  and  Thule,  Greenland),  where  on  some  occasions  the  predicted  values  were 
nearly  five  times  greater  than  the  measured  values.  On  the  other  hand,  there  was  mod¬ 
erately  good  agreement  between  the  measured  and  predicted  values  at  a  southern  high- 
latitude  site  (Arrival  Heights,  Antarctica).  The  best  agreement  was  observed  at  middle 
to  low  latitudes.  The  data  used  to  make  these  comparisons  were  measured  by  a  ELF/ 
VLF  measurement  system  (Fraser-Smith  et  al.,  1987). 

CCIR  Report  322-3  is  an  output  document  of  the  CCIR  XVIth  Plenary  Assembly 
held  in  Dubrovnik,  Yugoslavia,  in  1986.  It  was  produced  by  the  International  Working 
Party  (IWP)  6/2  of  CCIR  Study  Group  6.  It  was  based  on  ^e  work  of  Spaulding  and 
Washburn  (1985). 

Because  of  the  worldwide  acceptance  of  the  CCIR  Report  322-3  radio  noise  model, 
an  attempt  was  made  to  determine  the  cause  of  this  discrepancy.  It  was  found  that  the 
basis  for  this  error  was  in  the  procedure  used  to  prepare  the  measured  noise  data  for 
the  determination  of  a  global  numerical  representation  of  the  1-MHz  data.  Since  the 
contour  plots  of  the  1-MHz  radio  noise  in  CCIR  Report  322-3  were  in  turn  generated 
from  the  numerical  representation  thus  developed,  they  are  also  in  error. 

The  purpose  of  this  report  is  to  present  the  probable  cause  of  the  CCIR  322-3  ra¬ 
dio  noise  model.  First,  a  description  of  how  CCIR  Report  322  (CCIR,  1964)  noise 
models  have  been  and  are  currently  used  is  presented  to  provide  an  understanding  on 
how  a  systematic  error  in  a  1-MHz  noise  model  could  appear  at  other  frequencies.  Sec¬ 
ond,  the  development  process  of  the  current  CCIR  Report  322-3  model  is  presented  so 
that  the  cause  of  the  error  can  be  fully  understood.  Then  the  cause  itself  of  the  error 
in  the  1-MHz  model  along  with  the  frequency  dependence  of  this  error  for  three  meas¬ 
urement  sites  are  presented.  Next,  a  proposed  process  to  develop  a  corrected  1-MHz 
noise  model  is  discussed.  Finally,  conclusions  and  recommendations  are  made.  Primary 
among  these  recommendations  is  that  the  CCER  322-3  radio  noise  model  should  be 
used  with  caution  until  a  replacement  for  ^he  1-MHz  model  has  been  developed.  In 
particular  the  northern  and  southern  high  latitudes,  the  Arabian  Peninsula,  northern 
Africa,  and  the  mid-Atlantic  appear  to  be  the  most  inaccurate  areas.  Europe,  Asia,  the 
Indian  Ocesai,  Australia,  and  the  western  Pacific  from  Asia  to  the  date  line  appear  to 
be  the  most  accurate  areas. 
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2.0  THE  CCIR  322  NOISE  MODEL 

The  CCIR  Report  322  radio  noise  model  is  based  upon  measurements  of  average 
power  levels  (fa)  of  atmospheric  noise  on  a  worldwide  basis  starting  in  1957  by  the 
Central  Radio  Ftopagation  Laboratory  (CRPL).  Figure  1  shows  the  locations  of  the  16 
identical  recording  stations  used  to  measure  this  data.  At  most  of  the  sites  the  data 
continued  to  be  measured  until  1966.  Measurements  of  F,.  the  external  noise  figure 
defined  as  10  log  f*;  Vj.  the  dB  difference  between  the  average  voltage  and  the  rms 
voltage;  and  L^,  the  dB  difference  between  the  antilog  of  the  envelope  voltage  and  the 
rms  voltage,  were  made  using  a  bandwidth  of  200  Hz.  The  measurements  were  made 
on  frequencies  of  essentially  0.013,  0.051,  0.160,  0.495,  2.5,  5,  10,  and  20  MHz. 
Amplitude  probability  distribution  (APD)  measurements  were  made  at  some  of  the 
sites.  These  data  were  analyzed  and  published  in  a  series  of  reports  known  as  the  NBS 
Technical  Note  Series  18  (Crichlow  et  al.,  1959;  Crichlow,  Disney,  &  Jenkins,  1967). 
For  each  frequency  and  location,  the  month-hour  median  value  of  F*  along  with  Du 
and  Di,  the  upper  and  lower  decile  values,  were  given.  The  median  values  of  V^and 
Ld  were  also  given.  In  addition,  the  corresponding  season-time  block  values  of  these 
parameters  were  given  for  the  four  seasons,  winter  (December,  January,  and  Febru¬ 
ary),  spring  (March,  April,  May),  summer  (June,  July,  August),  and  fall  (September, 
October,  and  November),  and  six  4-hour  time  blocks  (0000-0400,  etc.).  The  atmos¬ 
pheric  noise  data  used  to  develop  CCIR  Report  322  were  the  data  from  July  1957 
through  October  1961. 
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Figure  1.  Radio  noise  recording  stations  used  to  obtain  data  used  to  develop 
the  original  CCIR  Report  322. 


The  CC3R  Report  322  radio  noise  model  consists  basically  of  two  parts:  a  1-MHz 
model  and  a  frequency-dependence  model.  Figure  2  shows  figure  5 A  of  CCIR  Report 
322.  This  figure  gives  Fam.  the  median  noise  figure,  at  1  MHz  as  a  function  of  latitude 
and  longitude  for  the  winter  season  and  the  time  block  1200-1600.  Since  this  map  is 
for  a  particular  season,  there  is  a  discontinuity  at  the  equator  (corresponding  to  winter 
being  6  months  apart  in  the  northern  and  southern  hemispheres).  To  obtain  Fan  at  a 
particular  frequency,  figure  3  (figure  SB  of  CCIR  Report  322)  is  used  to  convert  the  1- 
MHz  value  to  any  frequency  between  10  kHz  and  30  MHz.  Finally,  the  median  value 
of  Vd,  Vdmt  and  the  statistical  variations  of  Fa  about  its  median  value  Fam,  are  given 
also  in  figure  3  (figure  SC  in  CCIR  Report  322). 
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Figure  5A 

Expected  vatues  of  atmospheric  radio  noise.  r„„, 
(db  above  kTJb  at  1  Mc/s) 

(Winter:  1200-1600  h.) 


Figure  2.  Figure  5A  from  COR  Report  322. 
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Numerical  representation  of  CCIR  Report  322  is  available.  The  first  by  Lucas  and 
Harper  (1965)  covers  the  frequency  range  from  10  kHz  to  30  MHz.  This  representation 
was  obtained  by  numerically  mapping  values  obtained  from  the  CC3R  322  1-MHz  maps 
themselves,  rather  than  by  numerically  mapping  the  original  data  points  (84-longitude 
by  100-latitude  grid  points)  that  produced  the  CCIR  Report  322  maps.  Spaulding  and 
Washburn  (1985)  report  an  rms  error  of  2  dB  with  maximum  errors  up  to  approxi¬ 
mately  10  dB  for  the  Lucas  and  Harper  model  due  to  the  method  of  obtaining  the  rep¬ 
resentation.  The  Lucas  and  Harper  numerical  representation  of  the  frequency  variation 
of  F*ni,  Du,  and  Di  are  “precise,”  as  these  numerical  routines  were  used  to  produce 
these  particular  parts  of  CCIR  Report  322.  World  maps  of  1-MHz  atmospheric  noise  in 
universal  time  are  available  (Zacharisen  &  Jones,  1970).  These  universal  time  maps 
were  obtained  by  numerical  mapping  using  the  “original”  Report  322  data.  Also,  using 
the  original  data  to  plot  the  contour  maps  in  Report  322,  Sailors  and  Brown  (1982, 
1983)  developed  a  simplified  atmospheric  noise  model  suitable  for  use  on  minicom¬ 
puters. 
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3.0  THE  PROCESS  USED  TO  DEVELOP  THE  CCIR 
REPORT  322-3  RADIO  NOISE  MODEL 

CCIR  Report  322  was  developed  using  data  available  through  October  1961.  Since 
then  much  additional  data  has  become  available.  Data  from  the  original  worldwide  net¬ 
work  continued  to  be  measured  through  1966.  Many  years  of  data  from  10  Soviet 
measurement  locations  became  available  along  with  data  from  Thailand  for  March 
1966  through  February  1968  (Chindahpom  and  Younker,  1968).  All  of  this  data  was 
analyzed  and  updated  set  of  atmospheric  radio  noise  estimates  produced,  essentially  in 
the  CCIR  Report  322  format  (Spaulding  and  Washburn,  1985).  These  estimates  are 
new  1-MHz  contour  maps  corresponding  to  figure  2.  This  model  was  submitted  to  the 
CCIR  IWP  6/2  by  the  United  States  as  a  replacement  for  CCIR  Report  322-2  (CCIR 
Report  322  reprinted  with  a  revised  text  and  title,  but  with  the  same  atmospheric  noise 
estimates).  This  became  CCIR  Report  322-3. 

This  section  attempts  to  describe  the  process  used  to  develop  the  new  CCIR  Report 
322-3.  It  is  based  on  the  work  of  Spaulding  and  Washburn  (1985)  and  Spaulding 
(1992*). 

3.1  THE  NEW  DATA 

The  worldwide  network  locations  and  new  locations  are  given  in  table  1  (table  1  in 
Spaulding  and  Washburn  (1985)).  Figure  4  (figure  7  in  Spaulding  and  Washburn 
(1985))  is  a  repeat  of  figure  1,  but  with  the  new  locations  added. 

In  the  development  of  the  new  noise  model,  data  from  Thule,  Greenland,  and  Byrd 
Station,  Anarctica,  were  not  used.  These  data  were  not  used  because  it  was  assumed 
that  the  data  were  generally  contaminated  by  high  levels  of  man-made  noise. 

For  a  number  of  years,  the  Soviet  Union  operated  a  network  of  10  noise  measure¬ 
ment  stations.  Data  from  these  measurement  locations  within  the  old  Soviet  Union 
were  available  from  the  World  Data  Center  (National  Oceanic  and  Atmospheric  Ad¬ 
ministration,  Boulder,  CO).  Raw  data  were  available  on  microfilm  for  periods  of  time 
from  mid-1958  through  1965,  coincident  with  the  measurements  of  the  worldwide  net¬ 
work.  The  parameters  that  were  measured  were  different  from  those  discussed  for  the 
worldwide  network  above.  The  analysis  and  use  of  the  Soviet  data  are  discussed  next. 


Spaulding,  A.  D.,  private  communication,  May  1992. 


Table  1.  Atmospheric  noise  measurement  locations. 


WORLDWIDE  NETWORK  LOCATIONS  (CCIR  322^. 


Balboa,  Canal  Zone 

79.5W, 

9.0N 

Bill,  Wyoming 

105.2W, 

43.2N 

Boulder,  Colorado 

105.1W, 

40.  IN 

Byrd  Station,  Antarctia 

120.0W, 

80.0S 

Cook,  Australia 

130.4E, 

30.6S 

Enkoping,  Sweden 

17.3E, 

59.5N 

Front  Royal,  Virginia 

78.2W, 

38.8N 

Ibadan,  Nigeria 

3.9E, 

7.4N 

Kekaha,  Hawaii 

159.7W. 

22.0N 

New  Delhi,  India 

77. 3E, 

28.8N 

Ohira,  Japan 

140.5E, 

35. 6N 

Pretoria,  South  Africa 

28. 3E, 

25.8S 

Rabat,  Morocco 

6.8W. 

33.9N 

San  Jose,  Brazil 

45. 8W, 

23.3S 

Singapore 

103.8E, 

1.3N 

Thule,  Greenland 

68.7W, 

76.6N 

NEW  LOCATIONS 


Alma  Ata,  USSR 
Ashkhabad,  USSR 
Irkutsk,  USSR 
Khabarovsk,  USSR 
Kiev,  USSR 

Laem  Chabang,  Thailand 
Moscow,  USSR 
Murmansk,  USSR 
Simferopol,  USSR 
Sverdlovsk,  USSR 
Tbilisi,  USSR 


76.92E, 

43.25N 

8.3E, 

37.92N 

104.5E, 

52.0N 

135.0E, 

50.0N 

30. 3E, 

50.72N 

100.9E, 

13.05N 

37.32E, 

55.47N 

35.0e, 

69.0N 

34.03E, 

45.02N 

61.07E, 

56.73N 

40.0E, 

41.72N 

Figure  4.  Radio  noise  recording  locations  used  by  Spaulding  and  Washburn  (1985). 


3.2  ANALYSIS  OF  THE  SOVIET  DATA 

The  Soviet  atmospheric  noise  measurement  program  was  organized  and  controlled 
by  Dr.  Ja.  I.  Ukhter  (Izmiran,  P.  O.,  Akademgorodok,  Moscow  Region,  USSR).  Dr. 
Likhter  supplied  to  Dr.  A.  Donald  Spaulding  of  the  Institute  for  Telecommunication 
Sciences  in  Boulder,  Colorado,  detailed  information  on  the  measurement  equipment 
used  and  the  deflnitions  of  the  various  parameters  measured.  On  each  measurement 
frequency,  a  measurement  lasted  approximately  2  minutes,  and  measurements  were 
taken  3  hours  apart  each  day.  On  many  days,  no  measurements  were  made.  The 
voltage  levels  (given  in  field  strength,  pV/m)  that  were  exceeded  2,  10,  20,  30,  40,  50, 
60,  70,  80,  and  90  percent  of  the  time  were  recorded.  On  the  microfilm,  these  are 
noted  by  E0.02  to  E0.9.  Because  of  averaging  in  the  receiving  and,  perhaps,  the  short 
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measurement  time,  this  set  of  measurements,  unfortunately,  did  not  appear  to  Dr. 
Spaulding  to  be  APD  measurements,  and  they  did  not  seem  to  correspond  to  other  So¬ 
viet  APD  results.  Also,  most  of  the  energy  in  the  atmospheric  noise  process  is  con¬ 
tained  at  levels  that  occur  less  than  2  percent  of  the  time.  The  peak  value  was  also  re¬ 
corded  and  appears  on  the  microfilm.  Also  on  the  microfilm  is  another  parameter 
noted  as  Eon-  This  parameter  has  no  physical  meaning  in  itself  but  is  a  level  set  by 
the  equipment  operator,  below  which  the  other  levels  (E0.02.  etc.)  were  recorded. 
Spaulding  and  Washburn  (1985)  determined  that  Eon  would  serve  as  a  good  approxi¬ 
mation  to  the  rms  level  (fa),  and  the  parameter  Eon  was  used  by  them  in  the  analysis 
of  the  Soviet  data. 

The  median  value  of  Eon  was  determined  for  all  measurements  at  a  given  location 
and  for  a  given  measurement  frequency  for  the  hours  and  months  within  each  of  the 
twenty-four  3-month/4-hour  time  blocks.  This  median  value  of  Eon  is  in  pV/m.  The  an¬ 
tennas  used  at  the  Soviet  measurement  sites  were  5-meter  vertical  rods  over  a  ground 
plane.  To  go  from  field  strength  to  Fam,  therefore,  it  was  necessary  to  use  equation  1, 
the  equation  for  the  vertical  component  of  the  rms  field  strength  for  a  short  (h  «  X) 
grounded  vertical  monopole, 

En  =  F.  +  20  log  fMH2  +  B  -  95.5  dB  (1  pV/m)  (1) 

where  En  is  the  field  strength  [dB  (1  pV/m)]  in  bandwidth  b  (Hz),  f^Hz  is  the  center 
frequency  in  MHz,  and  B  =  10  log  b.  The  bandwidths  used  were  approximately  250  Hz 
at  frequencies  below  1.5  MHz,  and  1000  Hz  at  frequencies  at  or  above  1.5  MHz,  with 
some  variation  at  some  of  the  measurement  sites.  For  example,  a  1000-Hz  bandwidth 
was  sometimes  used  at  0.750  and  1  MHz.  Throughout  the  Soviet  data,  there  were 
missing  months,  time,  frequencies,  some  months  with  only  a  few  days  of  measure¬ 
ments.  All  in  all,  there  was  a  large  body  of  usable  data.  Some  of  the  data  was  ana¬ 
lyzed  by  The  Institute  for  Telecommunication  Sciences  (ITS)  in  Boulder,  Colorado,  but 
most  of  the  data  were  analyzed  by  the  Naval  Ocean  Systems  Center  (NOSC)*  in  San 
Diego,  California. 

The  measurement  frequencies  and  other  information  are  summarized  in  Spaulding 
and  Washburn  (1985)  for  each  of  the  Soviet  measurement  locations. 

The  analysis  involved  determining,  at  each  frequency,  for  each  3-month  period  and 
4-hour  time  block,  the  median  value  of  all  the  data.  These  median  values  at  the  vari¬ 
ous  frequencies  were  then  used  to  determine  the  approximate  1-MHz  Fan,  value.  This 
value  was  then  used  to  obtain  a  correction  value  to  the  CCIR  Report  322  value.  Figure 
5  (Egure  8  in  Spaulding  and  Washburn  (1985))  shows  an  example  for  Moscow  for 
June,  July,  August  (Northern  Hemisphere  summer)  and  1600-2000  hours.  Note  that  the 
noise  curve  for  a  quiet  receiving  location  falls  considerably  below  the  analyzed  median 
values.  A  computer  algorithm  was  developed  that  determined  the  atmospheric  noise 


*  Now  the  RDT&E  Division  of  the  Naval  Command,  Control  and  Ocean  Surveillance  Center. 
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frequency  variation  curve  that  “best”  fit  the  data.  However,  since  the  median  value  at 
some  frequencies  was  based  on  much  more  data  than  the  value  of  the  other  frequen¬ 
cies  for  a  location,  time  block,  and  season  (due  to  missing  data  and  some  frequencies 
being  stressed  at  some  locations);  “fitting”  process  was  generally  done  by  hand  (visu¬ 
ally).  On  figure  5,  the  “best”  fitting  frequency  law  curve  was  determined  to  be  72  dB. 
The  CCIR  Report  322  value  is  65  dB.  Hence,  a  value  of  +7  dB  can  be  used  to  correct 
the  atmospheric  noise  predicted  by  CCIR  Report  322  at  Moscow  during  the  summer  at 
1600-2000  hours.  Figure  6  (figure  9  in  Spaulding  and  Washburn  (1985))  shows  an  ex¬ 
ample  for  Moscow  for  the  period  November,  December,  January,  0800-1200  hours. 
Atmospheric  noise  would  be  expected  to  be  low  during  this  period  (winter  morning) 
and  could  possibly  be  contaminated  by  man-made  noise  at  the  higher  frequencies.  For 
the  higher  frequencies,  350  kHz  and  above,  the  figure  shows  a  typical  man-made  noise 
curve  at  a  level  to  be  expected  for  a  quiet  receiving  site.  Because  of  this  contamination 
possibility,  the  lower  frequencies  are  used  to  determine  the  frequency  law  curve.  In 
this  case,  the  frequency  law  curve  for  31  dB  was  determined.  The  CCIR  Report  322 
value  was  29  dB,  resulting  in  a  required  correction  of  +2  dB. 

3.3  CORRECTIONS  TO  CCIR  REPORT  322  1-MHz  Fam  VALUES 

The  procedure  illustrated  above  for  figures  5  and  6  for  determining  the  corrections 
to  be  made  to  CCIR  Report  322  was  followed  by  Spaulding  and  Washburn  (1985)  to 
obtain  corrections  for  each  location  and  for  each  3-month/4-hour  time  block.  Tables  2 
through  5  give  the  corrections  determined  using  this  procedure  for  each  of  the  four 
seasons.  Each  table  contains  a  correction  for  each  time  block  and  for  each  station  in 
table  1  except  for  those  noted  below.  The  “correction”  is  the  difference  between  the 
CCIR  Report  322  1-MHz  Fam  value  and  the  corresponding  value  determined  using  the 
above  procedure  from  the  data. 

For  certain  stations  listed  in  table  1,  no  corrections  were  determined.  No  correction 
values  were  obtained  for  Thule,  Greenland,  and  Byrd  Station,  Antarctica,  because  it 
was  assumed  that  the  data  were  contaminated  by  man-made  noise.  Since  there  were  no 
data  from  Ibadan,  Nigeria,  beyond  the  data  used  in  the  development  of  CCIR  Report 
322,  no  correction  was  used  for  Ibadan.  Because  the  corrections  for  Bill,  Wyoming, 
and  Boulder,  Colorado,  were  essentially  identical,  only  corrections  for  Boulder  were 
used.  Corrections  were  used  for  only  6  Soviet  locations  rather  than  10.  Simferopol, 
Sverdlovsk,  Tbilisi,  and  Kiev  either  had  only  small  amounts  of  usable  low-frequency 
data  necessary  to  determine  the  proper  1-MHz  F,m  value  or  were  close  to  other  meas¬ 
urement  locations.  The  data  at  these  four  Soviet  locations  were  analyzed  to  ascertain 
that  the  corrections  agreed  with  those  used  at  nearby  locations,  namely  Moscow  and 
Ashkhabad.  For  Murmansk  for  the  March,  April,  May  season  and  for  the  four  time 
blocks  0400-0800,  0800-1200,  1200-1600,  and  1600-2000  hours,  the  data  were  highly 
irregular  and  confusing.  Hence,  Spaulding  and  Washburn  (1985)  decided  not  to  at¬ 
tempt  to  obtain  any  correction  values  for  Murmansk  for  these  four  periods  (see 
table  3). 
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(dB  above  kTob) 


Figure  5.  Determination  of  1-MHz  Fan»  value  for  Moscow,  June, 
July,  August,  1600-2000  hours  (Spaulding  and  Washburn,  1985). 
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Figure  6.  Determination  of  1-MHz  F»ni  value  for  Moscow,  December, 
January,  February,  1000-1200  hours  (Spaulding  and  Washburn,  1985). 


Because  the  original  contour  maps  in  CQR  Report  322  were  produced  directly  from 
a  grid  of  equally  spaced  84-Iongitude  by  100-latitude  points  for  each  time  block/season, 
the  next  step  in  the  analysis  was  to  translate  the  correction  data  for  the  19  sites  to  this 
same  84-by-lOO  lattice.  Then  data  for  a  new  noise  model  were  found  simply  by  adding 
the  corrections  point-by-point  for  the  grid  point  of  original  data. 

To  do  this,  Spaulding  and  Washburn  (1985)  used  an  interpolation  method  due  to 
Dr.  Charles  L.  Lawson  (1982,  1984)  for  interpolating  scattered  data  over  a  sphere. 

This  method  first  constructs  a  triangular  grid  over  a  surface  (in  this  case  the  Earth) 
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using  a  given  set  of  points  as  vertices  (in  this  case  the  19  data  points  in  tables  2 
through  5).  Then  continuous  first  partial  derivatives  are  estimated  by  the  method  at 
each  vertex  using  local  quadratic  least  squares  approximations  to  given  data  values  at 
nearby  vertices.  The  method  for  interpolation  then  uses  six  Hermite  cubic  interpola¬ 
tions  along  arcs  of  great  circles.  Figures  7  through  30  are  the  resultant  contour  maps 
of  the  24  (four  3-month  periods,  six  4-hour  time  blocks)  lOO-by-84  correction  grids 
produced  by  Spaulding  and  Washburn.  These  maps  show  the  changes  to  the  CCIR 
Report  322  to  be  made.  There  are  substantial  corrections  in  some  areas  as  could  also 
be  seen  by  examining  tables  2  through  5.  The  correction  maps  are  presented  in  terms 
of  3-month  periods  rather  than  as  season  (which  results  in  a  discontinuity  at  the  equa¬ 
tor)  as  in  CCIR  Report  322. 


Table  2.  Corrections  (dB)  to  CCIR  Report  322  1-MHz  Fjun  values  for 
December,  January,  and  February  (Spaulding  and  Washburn  (1985)). 


PLACE 

LOCATION 

00-04 

04-08 

LOCAL  TIME 

08-12  12-16 

16-20 

20-24 

Alma  Ata 

76.9E,  43.2N 

-7 

-6 

+6 

+5 

-3 

-6 

Irkutsk 

104.5E,  52.0N 

-21 

-25 

-7 

-15 

-25 

-25 

Khabarovsk 

135.0E,  50.0N 

-19 

-15 

-8 

-7 

-20 

-20 

New  Delhi 

77.3E,  28.8N 

-13 

+7 

+17 

+17 

+8 

+11 

Ohira 

140.5E,  35.6N 

+7 

+7 

+8 

+12 

+11 

11 

Thailand 

100.9E,  13.0N 

+14 

+15 

+24 

+18 

+17 

+15 

Singapore 

103.8E,  1.3N 

0 

+6 

12 

+9 

+5 

+1 

Kekaha 

159.7W,  22.0N 

+5 

+10 

+8 

+15 

+5 

+5 

Boulder 

105.1W,  40.1N 

+5 

+4 

+7 

+14 

+7 

+8 

Front  Royal 

78.2W,  38.8N 

-1 

+2 

+3 

+8 

0 

0 

Balboa 

79.5W,  9.0N 

+3 

+6 

+7 

+9 

+7 

+2 

Rabat 

6.8W,  33.9N 

+2 

+4 

+3 

+8 

+2 

+4 

&ikoping 

17.3E,  59.5N 

+12 

+10 

-1 

+8 

+7 

+7 

Murmansk 

35.0E,  69.0N 

+8 

+5 

+7 

+9 

+7 

+7 

Moscow 

37.3E,  55. 5N 

+4 

+3 

+2 

+4 

0 

-1 

Ashkabad 

58.3E,  37.9N 

-9 

-1 

-5 

-5 

-6 

-12 

Cook 

130.4E,  30.6S 

+2 

-3 

+6 

+1 

+6 

+3 

San  Jose 

45. 8W,  23. 3S 

+2 

0 

+2 

+2 

+4 

+3 

Piretoria 

28.3E,  25.8S 

-4 

+8 

-4 

+1 

+5 

-8 
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Table  3.  Corrections  (dB)  to  CCIR  Report  322  l-MHzFam  values  for  March, 
April,  and  May  (Spaulding  and  Washburn  (1985)). 


PLACE 

LOCAL  TIME 

LOCATION 

00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Alma  Ata 

76.9E,  43.2N 

-7 

-4 

-5 

-8 

-2 

-5 

Irkutsk 

104.5E,  52.0N 

-12 

-7 

-5 

+5 

+4 

-14 

Khabarovsk 

135.0E,  50.0N 

-15 

-6 

-7 

+1 

-1 

-21 

New  Delhi 

77.3E,  28.8N 

+6 

+10 

+15 

+9 

+12 

+7 

Ohira 

140.5E,  35.6N 

+2 

+4 

+15 

+12 

+7 

+2 

Thailand 

100.9E,  13.0N 

+6 

+9 

+14 

+17 

+10 

+8 

Singapore 

103.8E,  1.3N 

+3 

+5 

+16 

+13 

+10 

+5 

Kekaha 

159.7W,  22.0N 

+6 

+8 

+11 

+13 

+5 

+6 

Boulder 

105.  IW,  40.  IN 

+3 

+5 

+7 

+2 

+8 

+3 

Front  Royal 

78.2W.  38.8N 

-3 

-1 

-5 

-5 

-3 

+1 

Balboa 

79.5W,  9.0N 

+4 

+5 

+9 

+5 

+6 

+4 

Rabat 

6.8W,  33.9N 

+1 

+4 

+3 

-5 

0 

0 

Enkoping 

17.3E,  59.5N 

0 

0 

+3 

-1 

-5 

+2 

Murmansk 

35.0E,  69.0N 

+3 

+4 

Moscow 

37.3E,  55.5N 

+4 

0 

0 

-2 

0 

-4 

Ashkabad 

58.3E,  37.9N 

+2 

+1 

-3 

+2 

-4 

+2 

Cook 

130.4E,  30.6S 

+3 

+9 

+5 

+8 

+4 

+5 

San  Jose 

45.8W.23.3S 

+2 

+3 

+5 

+7 

+3 

+7 

Pretoria 

28.3E,  25.8S 

+3 

+2 

+11 

+9 

+11 

-1 
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Table  4.  Corrections  (dB)  to  CCIR  Report  322  1-MHz  Fam  values  for  June, 
July,  and  August  (Spaulding  and  Washburn  (1985)). 


PLACE  LOCATION 


LOCAL  TIME 


00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Alma  Ata 

76.9E,  43.2N 

-4 

0 

-8 

-6 

-2 

-3 

Irkutsk 

104.5E,  52.0N 

-20 

-6 

-11 

0 

-4 

-15 

Khabarovsk 

135.0E,  50.0N 

-10 

-4 

-8 

+1 

+2 

-12 

New  Delhi 

77.3E,  28.8N 

+8 

+17 

+11 

+4 

+10 

+8 

Ohira 

140.5E,  35.6N 

+2 

+5 

+11 

+10 

+9 

+3 

Thailand 

100.9E,  13.0N 

+11 

+15 

+15 

+18 

+13 

+8 

Singapore 

103.8E,  1.3N 

+4 

+11 

+15 

+15 

+10 

+2 

Kekaha 

159.7W,  22.0N 

+9 

+6 

+2 

+2 

-6 

+4 

Boulder 

105.1W,40.1N 

+2 

+8 

+7 

+10 

+12 

+6 

Front  Royal 

78.2W,  38.8N 

-8 

-4 

+4 

-11 

-10 

-1 

Balboa 

79.5W.  9.0N 

-10 

+9 

+12 

+1 

+3 

+4 

Rabat 

6.8W,  33.9N 

0 

+3 

+2 

-16 

-4 

0 

Enkoping 

17.3E,  59.5N 

+5 

+6 

-4 

-7 

-4 

-7 

Murmansk 

35.0E,  69.0N 

-2 

+8 

-1 

+5 

+10 

-2 

Moscow 

37.3E,  55.5N 

-2 

0 

-2 

+2 

+7 

-6 

Ashkabad 

58.3E,  37.9N 

+2 

-4 

-7 

-8 

-3 

-3 

Cook 

130.4E,  30.6S 

+5 

+7 

+12 

+11 

+10 

+6 

San  Jose 

45.8W,  23.3S 

-4 

+5 

+11 

+10 

+9 

0 

Pretoria 

28.3E,  25.8S 

+12 

+11 

+20 

+22 

+16 

+17 
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Table  5.  Corrections  (dB)  to  CCIR  Report  322  1-MHz  Fam  values  for 
September,  October,  and  November  (Spaulding  and  Washburn  (1985)). 


PLACE 

LOCATION 

00-04 

Alma  Ata 

76.9E,  43.2N 

-4 

Irkutsk 

104.5E,  52.0N 

-22 

Khabarovsk 

135.0E,  50.0N 

-19 

New  Delhi 

77.3E,  28.8N 

+5 

Ohira 

140.5E,  35. 6N 

+6 

Thailand 

100.9E,  13.0N 

+5 

Singapore 

103.8E,  1.3N 

+7 

Kekaha 

159.7W,  22.0N 

+1 

Boulder 

105.1W,  40.1N 

+2 

Front  Royal 

78.2W,  38.8N 

-2 

Balboa 

79.5W,  9.0N 

0 

Rabat 

6.8W,  33.9N 

+5 

Enkoping 

17.3E,  59.5N 

+2 

Murmansk 

35. OE,  69.0N 

-5 

Moscow 

37.3E,  55.5N 

0 

Ashkabad 

58.3E.  37.9N 

+3 

Cook 

130.4E,  30.6S 

-1 

San  Jose 

45.8W,23.3S 

+4 

Pretoria 

28.3E,  25.8S 

+3 

LOCAL  TIME 


04-08 

08-12 

12-16 

16-20 

20-24 

-3 

-2 

-8 

-7 

-9 

-20 

-15 

-15 

-20 

-20 

-10 

-8 

-9 

-12 

-18 

•1-8 

+9 

-4 

+6 

+5 

+4 

+12 

+9 

+9 

+7 

+11 

+20 

+12 

+9 

+7 

+11 

+20 

+14 

+7 

+7 

+5 

0 

+2 

0 

+2 

+7 

+12 

+10 

+9 

+3 

+3 

+4 

-1 

-2 

-1 

+4 

+14 

+12 

+5 

0 

+9 

+11 

+3 

+6 

+5 

+4 

0 

+4 

0 

+3 

+2 

-2 

+3 

+1 

-2 

+3 

+3 

-2 

-2 

+2 

+5 

-6 

-4 

-2 

0 

+2 

+11 

+10 

+4 

+4 

+6 

+12 

+6 

+3 

+2 

+9 

+9 

+6 

+8 

+4 
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SMlBggaW—MMWBBMWiag 


75  00  105  120  135  150  155  ISO  105  150  135  120  105  90  75  60  45  30  15  0  15  ^  «5  > 

Figure  7.  Corrections  (dB)  to  original  CQR  Report  322  1-MHz  F*ni  estimates, 
December,  January,  February,  0000-0400  hours  (Spaulding  and  Washburn,  1985) 


5  90  ’05  17?  135  '5C  105  19C  10S  150  135  120  105  9C  75  60  45  30  15  0  15  30  45  00 
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Figure  10.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  estimates, 
Dumber,  January,  February,  1200-1600  hours  (Spaulding  and  Washburn,  1985). 
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Figure  11.  Corrections  (dB)  to  original  CCIR  Report  322  l>MHz  Fun  estimates, 
December,  January,  February,  1600-2000  hours  (Spaulding  and  Washburn,  1985) 
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Figure  12.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F^n  estimates. 
Dumber,  January,  February,  2000-2400  hours  (Spaulding  and  Washburn,  1985) 
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Figure  14.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F«ni  estimates, 
March,  April,  May,  0400-0800  hours  (Spaulding  and  Washburn,  1985). 
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•0  79  60  109  120  139  190  169  160  169  190  139  120  1C»  60  79  60  49  30  15  0  19  30  49  60 

Figure  IS.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F«n>  estimates, 
March,  April,  May,  0800-1200  hours  (Spaulding  and  Washburn,  1985). 
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Hgure  16.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  ^•in  CStllYlStCSf 
Mirch,  April,  May,  1200-1600  hours  (Spaulding  and  Washburn,  1985). 
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Figure  17.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  Fam  estimates, 
March,  April,  May,  1600-2000  hours  (Spaulding  and  Washburn,  1985). 


figure  18.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  Fam  estimates, 
N^ch,  April,  May,  2000-2400  hours  (Spaulding  and  Washburn,  1985). 
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Figure  19.  Corrections  (dB)  to  original  CCIR  Report  322  l-MHz  estimates, 
June,  July,  August,  0000-0400  hours  (Spaulding  and  Washburn,  1985). 


Figure  20.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F^m  estimates, 
June,  July,  August,  0400-0800  hours  (Spaulding  and  Washburn,  1985). 
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Figure  21.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F^nj  estimates, 
June,  July,  August,  0800-1200  hours  (Spaulding  and  Washburn,  1985). 
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Figure  22.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  Fun  estimates, 
June,  July,  August,  1200-1600  hours  (Spaulding  and  Washburn,  1985). 
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Figure  23.  Corrections  (dB)  to  original  COR  Report  322  1-MHz  Fua  estimates, 
June,  July,  August,  1600-2000  hours  (Spaulding  and  Washburn,  198S). 


Figure  24.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  Fam  estimates, 
June,  July,  August,  2000-2400  hours  (Spaulding  and  Washburn,  1985). 
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Figure  25.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  estimates, 
September,  October,  November,  0000-0400  hours  (Spaulding  and  Washburn,  1985) 
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Figure  26.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F^n  estimates, 
September,  October,  November,  0400-0800  hours  (Spaulding  and  Washburn,  1985). 


Figure  27.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F,n»  estimates, 
September,  October,  November,  0800-1200  hours  (Spaulding  and  Washburn,  1985). 
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Hfure  28.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F«ni  estimates, 
September,  October,  November.  1200-1600  hours  (Spaulding  and  Washburn,  1985). 
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Figure  29.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  estimates, 
September,  October,  November,  1600-2000  hours  (Spaulding  and  Washburn,  1985) 
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Figure  30.  Corrections  (dB)  to  original  CCIR  Report  322  1-MHz  F.n,  estimates, 
September.  October,  November,  2000-2400  hours  (Spaulding  and  Washburn,  1985). 
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3.4  THE  NEW  l-MHz  Fam  VALUES 


The  new  1-MHz  Fam  data  values  for  constructing  a  new  1-MHz  noise  model  were 
obtained  by  Spaulding  and  Washburn  (1985)  by  adding  each  of  the  84-longitude  by 
100-latitude  grid  of  correction  values  to  the  corresponding  original  84-longitude  by  100- 
latitude  grid  data  values  from  which  CCIR  Report  322  was  constructed. 

In  developing  a  numerical  representation  for  these  new  1-MHz  Fa®  maps,  Spauld¬ 
ing  and  Washburn  used  the  method  used  by  Lucas  and  Harper  (1965).  The  resulting 
sets  of  numerical  coefficients  could  then  be  used  with  existing  programs  developed  to 
obtain  the  1-MHz  Fam  noise  value  from  the  Lucas  and  Harper  representation  of  the 
CCIR  Report  322  noise  model.  Spaulding  and  Washburn  give  details  on  how  these  co¬ 
efficients  were  obtained. 

Spaulding  and  Washburn  (1985)  compared  the  numerical  representation  thus 
obtained  for  each  of  the  8400  original  data  points  (84  x  100  grids)  for  each  of  the  24 
numerical  maps.  They  found  an  rms  variation  that  ranged  from  0.88  dB  to  2.37  dB 
over  the  24  maps  with  an  average  rms  variation  of  1.52  dB  and  with  a  maximum  de¬ 
viation  of  6.7  ®  (all  maps  considered,  i.e.,  24  x  8400  points). 

These  numerical  maps  represent  a  “smoothed”  version  of  the  original  data  and  are 
the  new  1-MHz  F^,  worldwide  atmospheric  noise  estimates.  Figures  31  through  54  (fig¬ 
ures  34  through  57  from  Spaulding  and  Washburn  (1985))  are  contour  plots  of  these 
estimates. 

The  contour  plots  in  CCIR  Report  322-3  are  similar,  except  the  CCIR  I\VT  6/2 
gathered  the  plots  together  by  season  rather  than  by  months  as  in  figures  31  through 
54.  This  resulted  in  the  discontinuity  at  the  equator  in  CCIR  Report  322-3. 
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Figure  31.  Expected  values  of  atmospheric  radio  noise  at  1  MHz, 

(dB  above  kTob),  for  December,  January,  February,  0000-0400  hours 
(Spaulding  and  Washburn,  1985). 
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Hgure  32.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F^ 
(<ffl  above  kTob),  for  December,  January,  February,  0400-0800  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  33.  Expected  values  of  atmospheric  radio  noise  at  1  MHz, 

(dB  above  kTob),  for  December,  January,  February,  0800-1200  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  34.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F»n, 
(dB  above  kTob),  for  December,  January,  February,  1200-1600  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  35.  Expected  values  of  atmospheric  radio  noise  at  1  MHz, 

(dB  above  kTob),  for  December,  January,  February,  1600-2000  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  36.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fam 
(dB  above  kTob),  for  December,  January,  February,  2000-2400  hours 
(Spaulding  and  Washburn,  1985). 


33 


90  -05  rrc  1j''  'Sr  'as  tso  fas  iSO  I35  »20  ’O*  ©'•  ’5  so  30 


jRfl8fri8iSIS08^^^S98 
■^^ar/iRRSiaMfSTi 


iptjii 


79  to  10s  120  t3S  ISO  las  180  t<s  ISO  13S  120  10S 


Figure  37.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fu 
(dB  above  kTob),  for  March,  April,  May,  0000-0400  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  38.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fu, 
(tffi  above  kTob),  for  March,  April,  May,  0400-0800  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  39.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F^n 
(dB  above  kTob).  for  March,  April,  May,  0800-1200  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  40.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fu 
(dB  above  kTob),  for  March,  April,  May,  1200-1600  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  41.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fu, 
(dB  above  kTob),  for  March,  April,  May,  1600-2000  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  42.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fjn 
(<S  above  kTob),  for  March,  April,  May,  2000-2400  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  44.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F„ 
(dB  above  kTob),  for  June,  July,  August,  0400-0800  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  45.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F»ni 
(dB  above  kTob),  for  June,  July,  August.  0800-1200  hours  (Spaulding 
and  Washburn,  1985). 


Figure  46.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F»n> 
(<ffl  above  kT^),  for  June,  July,  August,  1200-1600  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  47.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F*, 
(dB  above  kTob),  for  June,  July,  August,  1600-2000  hours  (Spaulding 
and  Washburn,  1985). 
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Figure  48.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F„ 
(dB  above  kTob),  for  June,  July,  August,  2000-2400  hours  (Spaulding 
and  Washburn,  1985). 


Figure  49.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F^n 
(dB  above  kTob),  for  September,  October,  November,  0000-0400  hours 
(Spaulding  and  Washburn,  1985). 
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Hgure  50.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F^nj 
(dB  above  kTob),  for  September,  October,  November,  0400-0800  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  51.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  F^n 
(dB  above  kTob),  for  September,  October,  November,  0800-1200  hours 
(Spaulding  and  Washburn,  1985). 
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Figure  52.  Expected  values  of  atmospheric  radio  noise  at  1  MHz,  Fun 
(<ffl  above  kTob),  for  September,  October,  November,  1200-1600  hours 
(Spaulding  and  Washburn,  1985). 
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4.0  THE  PROBABLE  SOURCE  OF  THE  ERROR  IN  THE 
CCIR  RADIO  NOISE  MODEL 


The  probable  source  of  the  error  in  the  CQR  radio  noise  model  is  most  likely  due 
to  the  nonuse  of  correction  factors  for  Thule,  Greenland;  Byrd  Station,  Antarctica; 
Ibadan,  Nigeria;  and  Bill,  Wyoming.  In  the  case  of  Thule  and  Byrd  Station,  data  from 
these  sites  were  not  used  because  of  possible  contamination  from  man-made  noise  at 
these  two  locations.  However,  these  same  two  locations  were  used  in  the  original  CGR 
noise  model.  One  would  anticipate  that  if  the  data  were  contaminated  by  man-made 
noise  that  a  negative  correction  factor  would  be  the  likely  case.  Thus,  it  would  seem 
desirable  to  determine  correction  factors  for  these  sites  as  was  done  for  the  other  sites. 
That  is,  the  data  at  lower  frequencies  could  have  been  used  to  obtain  correction  factors 
at  1  MHz.  Certainly,  there  is  evidence  that  the  measurements  at  Thule  at  2.5  and  5.0 
MHz  were  contaminated  by  man-made  noise  (Herman,  1962).  However,  Herman 
(1963)  showed  that  Byrd  Station  is  an  exceptionally  quiet  location  and  would  be  a 
good  site  for  making  a  variety  of  radio  measurements  requiring  a  low-noise  back¬ 
ground.  At  Thule,  man-made  noise  on  2.5  and  5.0  MHz  appeared  to  be  57  and  49  dB 
above  kTB,  respectively,  while  at  Byrd  Station  the  values  were  about  20  and  12  dB, 
respectively  (Herman,  1964).  These  values  were  estimated  from  data  taken  during  a 
PCA,  when  atmospheric  noise  was  absent.  The  corresponding  values  for  a  quiet  rural 
site  at  2.5  and  5.0  MHz  are  43  and  34  dB,  respectively.  Both  sites  are  affected  by  ga¬ 
lactic  noise  at  10.0  and  20.0  MHz.  In  the  case  of  Ibadan,  no  additional  data  were 
available  beyond  the  last  date  for  data  used  in  the  original  CCIR  322  noise  model.  Not 
even  a  correction  factor  of  zero  to  maintain  the  status  quo  was  used  for  these  three 
sites.  In  the  case  of  Bill,  Wyoming,  data  were  not  included  because  the  correction  fac¬ 
tors  obtained  for  it  were  nearly  the  same  as  for  Boulder,  Colorado,  which  was  in  close 
proximity  to  Bill,  Wyoming. 


The  net  result  of  not  including  correction  factors  (not  even  zero  to  maintain  the 
status  quo)  for  these  four  locations  was  that  the  interpolation  algorithm  used  to  deter¬ 
mine  the  100-latitude  by  84-Iongitude  data  points  supplied  erroneous  correction  factors 
at  these  sites.  Tables  6  through  9  give  the  error  in  the  correction  factor  for  these  four 
sites.  For  Thule,  Byrd  Station,  and  Ibadan,  the  error  is  the  difference  between  the  cor¬ 
rection  factors  given  in  figures  7  through  30  and  zero  for  the  status  quo.  For  Bill, 
Wyoming,  the  error  is  the  difference  between  those  given  in  these  figures  and  the  cor¬ 
rection  factor  input  for  Boulder,  Colorado.  It  was  anticipated  that  the  errors  at  Thule 
and  Byrd  Station  might  be  large,  but  it  was  a  surprise  to  see  the  magnitude  of  the  er¬ 
rors  at  Ibadan.  For  Thule  the  maximum  and  minimum  errors  in  the  correction  con¬ 
tours  from  zero  status  quo  were  10.1  and  -10.8  dB,  respectively.  For  Ibadan  the  maxi¬ 
mum  and  minimum  errors  were  12.5  and  -1.5  dB,  respectively.  For  Byrd  Station  the 
maximum  and  minimum  errors  were  12.0  and  3  dB,  respectively.  The  error  for  Bill, 
Wyoming,  is  within  the  rms  error  of  the  numerical  maps  of  CCIR  Report  322-3. 
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Table  6.  Interpolation  errors  (dB)  for  select  measurement  locations 
for  December,  January,  and  February. 


LOCATION  NAME  LOCAL  TIME 


00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Thule 

7.0 

3.5 

6.0 

6.0 

3.0 

6.0 

Ibadan 

-0.4 

6.0 

0.6 

5.5 

3.0 

-1.5 

Bill 

-0.8 

-0.6 

-1.0 

-0.3 

-1.0 

-1.0 

Byrd  Station 

3.0 

0.0 

5.5 

2.0 

6.0 

5.5 

Table  7.  Interpolation 

errors  (dB)  for  select  measurement 

locations 

for  March,  April,  and 

May. 

LOCATION  NAME 

LOCAL  TIME 

00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Thule 

3.0 

-3.0 

-3.0 

-4.1 

-10.8 

3.0 

Ibadan 

3.0 

3.0 

8.2 

0.8 

5.0 

0.0 

Bill 

-1.0 

-0.9 

-1.0 

-0.5 

-1.0 

-0.6 

Byrd  Station 

4.0 

8.0 

12.0 

7.3 

10.2 

7.0 

Table  8.  Interpolation 

errors  (dB)  for  select  measurement 

locations 

for  June,  July,  and  August. 

LOCATION  NAME 

LOCAL  TIME 

00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Thule 

-1.3 

6.0 

1.5 

6.0 

10.1 

0.0 

Ibadan 

5.3 

7.4 

10.6 

0.0 

4.5 

8.0 

Bill 

-0.5 

-1.0 

-1.0 

-1.0 

-1.0 

-0.7 

Byrd  Station 

8.3 

8.0 

12.0 

7.3 

10.2 

7.0 

Table  9.  Interpolation  errors  (dB)  for  select  measurement  locations 
for  September,  October,  November. 


LOCATION  NAME  LOCAL  TIME 


00-04 

04-08 

08-12 

12-16 

16-20 

20-24 

Thule 

-6.0 

3.0 

0.0 

3.0 

3.0 

-2.5 

Ibadan 

5.5 

10.7 

12.5 

5.2 

8.8 

5.3 

Bill 

-1.4 

-1.0 

-1.0 

-1.0 

-0.7 

-0.6 

Byrd  Station 

3.0 

4.9 

12.0 

9.0 

3.0 

3.0 
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4.1  GEOGRAPHICAL  EFFECTS 

To  see  the  geographical  effects  of  these  errors,  contour  plots  were  made  of  the  er¬ 
rors  for  each  time  block  and  season.  For  the  19  locations  used  to  determine  the  correc¬ 
tion  factor  contours  in  figures  7  through  30,  the  error  was  assumed  to  be  zero.  That  is, 
it  was  assumed  that  the  interpolation  process  gives  values  for  the  correction  factors  at 
these  sites  as  were  input  into  the  interpolation  process.  The  values  given  in  tables  6 
through  9  were  used  for  the  other  four  sites.  The  commercial  graphics  program  for 
personal  computers,  Axum\  was  used  to  plot  the  errors  for  the  irregularly  spaced  data 
points.  Axum  rarely  was  able  to  determine  a  grid  spacing  greater  than  50  by  50  for  its 
internal  interpolation  for  the  contour  plotting.  The  results  are  presented  in  figures  55 
through  78.  Note  that  the  longitudes  are  positive  for  degrees  east  of  Greenwich  (zero 
degree)  and  negative  for  degrees  west  of  zero  degree.  This  longitude  convention  was 
used  so  that  the  contour  plots  produced  by  Axum  would  have  the  eastward  direction 
on  the  right-hand  side  as  is  the  convention  for  maps.  Postive  latitudes  are  northern  lati¬ 
tudes,  and  negative  latitudes  are  southern  latitudes.  The  locations  of  the  four  sites  for 
which  no  correction  factors  were  used  in  the  interpolation  are  approximately  given  in 
the  flgures.  Examination  of  these  figures  reveals  ^at  the  geographical  extent  of  the 
error  is  not  confined  to  the  measurement  location  but  in  fact  in  some  cases  is  very 
large.  This  is  especially  true  in  the  northern  and  southern  high  latitudes,  the  Arabian 
Peninsula,  northern  Africa,  and  the  mid-Atlantic  areas. 

4.2  FREQUENCY-DEPENDENCE  EFFECTS 

As  the  CCIR  Report  322-3  atmospheric  noise  model  can  be  used  at  frequencies 
from  10  kHz  through  30  MHz,  it  is  important  to  be  able  to  translate  an  error  in  the  1- 
MHz  model  to  an  error  at  any  frequency  in  this  frequency  range.  To  determine  this 
relationship,  the  frequency-dependence  model  was  used,  assuming  no  error  inherit  in 
it.  This  model  obtains  the  noise  at  any  frequency  by  inputing  the  1-MHz  value  for  a 
particular  time  block  and  season  into  a  set  of  curves  parametric  in  the  1-MHz  value 
(see  figure  3),  For  each  time  block  and  season,  there  is  a  range  of  possible  input  1- 
MHz  values  for  which  the  atmospheric  noise  can  be  obtained  depending  on  the  loca¬ 
tion  of  the  receive  site.  To  obtain  the  error  at  an  arbitrary  frequency,  time  block,  and 
season,  each  of  the  parametric  values  was  perturbed  by  errors  of  10  and  -10  dB,  and 
the  resultant  errors,  respectively,  were  determined.  The  statistical  average,  high,  and 
low  values  were  then  obtained.  It  was  found  that  the  error  at  each  of  the  35  frequen¬ 
cies  from  10  kHz  through  30  MHz  used  did  not  depend  significantly  on  the  parametric 
curve  value.  However,  the  error  was  diumally  dependent.  The  maximum  errors  for 
each  frequency  occurred  during  the  local  time  daytime,  and  the  minimum-error  errors 
occurred  during  the  nighttime.  Figure  79  shows  the  average,  high,  and  low  error  as  a 
function  of  frequency  and  1-MHz  error.  The  figure  shows  the  error  for  -10  dB  to  be  a 
mirror  image  of  that  for  a  10-dB  error. 


^  trademark  of  TrlMetrix,  Inc. 
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Figure  55.  Geographical  l-MHz  atmospheric  noise  model  error,  December,  January, 
February,  0000-0400  hours. 


Figure  56.  Geographical  1-MHz  atmospheric  noise  model  error,  December,  January, 
February,  0400-0800  hours. 
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Figure  57.  Geographical  1-MHz  atmospheric  noise  model  error,  December, 
January,  February,  0800-1200  hours. 


Figure  58.  Geographical  1-MHz  atmospheric  noise  model  error,  December, 
January,  February,  1200-1600  hours. 
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Figure  59.  Geographical  l-MHz  atmospheric  noise  model  error,  December, 
January,  February,  1600-2000  hours. 
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Figure  60.  Geographical  1-MHz  atmospheric  noise  model  error,  December, 
January,  February,  2000-2400  hours. 
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Figure  61.  Geographical  l-MHz  atmospheric  noise  model  error,  March, 
April,  May,  0000-0400  hours. 


Figure  62.  Geographical  1-MHz  atmospheric  noise  model  error,  March, 
^ril,  May,  0400-0800  hours. 
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Figure  65.  Geographical  1-MHz  atmospheric  noise  model  error,  March, 
April,  May,  1600-2000  hours. 
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Figure  66.  Geographical  1-MHz  atmospheric  noise  model  error,  March, 
^ril.  May,  2000-2400  hours. 
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Figure  69.  Geographical  1-MHz  atmospheric  noise  model  error,  June, 
July,  August,  0800-1200  hours. 
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Figure  70.  Geographical  l>MHz  atmospheric  noise  model  error,  June, 
July,  August,  120(1^1600  hours. 
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Figure  72.  Geographical  1-MHz  atmospheric  noise  model  error,  June, 
July,  August,  2000-2400  hours. 
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Figure  73.  Geographical  1-MHz  atmospheric  noise  model  error,  September, 
O^ber,  November,  0000-0400  hours. 


Figure  74.  Geographical  l-MHz  atmospheric  noise  model  error,  September, 
October,  November,  0400-0800  hours. 
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Figure  75.  Geographical  1-MHz  atmospheric  noise  model  error,  September, 
October,  November,  0800-1200  hours. 


Figure  76.  Geographical  1-MHz  atmospheric  noise  model  error,  September, 
October,  November,  1200-1600  hours. 
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Figure  77.  Geographical  1-MHz  atmospheric  noise  model  error,  September, 
October,  November,  1600-2000  hours. 
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Figure  78.  Geographical  1-MHz  atmospheric  noise  model  error,  September, 
O^ber,  November,  2000-2400  hours. 
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Figure  79.  Error  in  atmospheric  noise  as  a  function  of  frequency  for  a  given 
error  at  1  MHz,  mean,  maximum,  minimum  for  all  time  blocks. 


4.3  INTERPOLATION  EFFECTS 

The  accuracy  of  the  interpolation  itself  is  affected  by  the  lack  of  inclusion  of  data 
from  the  four  measurement  locations.  First,  the  number  of  triangles  and  number  of 
edges  of  these  triangles  is  reduced.  From  Lawson  (1977),  the  number  of  triangles  is 

ni  «  nb  +  2  (ni  -  1)  S  2n  (2) 

and  the  number  of  edges  is 

n,  «  2nb  +  3  (ni  -  1)  S  3n  (3) 

where  n  is  the  number  of  distinct  points  in  the  set  S;  nb  denotes  the  number  of  points 
in  S  on  the  boundary  of  the  convex  hull  of  S;  and  nj  denotes  the  number  of  points  in 
the  interior  of  the  convex  hull  of  S  so  that  n  «  nb  ■«'  ni.  In  the  case  of  data  points  on  a 
sphere,  nb  ■  0  (i.e.,  there  are  no  boundary  points).  For  n  »  19,  nt  *  36  and  n,  =  54.  For 
n  >  23,  nt  ■  44  and  n«  >  66.  Thus  there  is  a  122.22-percent  increase  in  triangle  edges 
by  adding  the  four  points.  Second,  the  number  of  input  data  points  affects  how  the 
surface  interpolation  determines  the  triangles  themselves.  Each  pair  of  triangles  forms 
a  quadrilateral.  The  Lawson  criterion  gives  the  preferred  triangulation  of  a  quadrilat¬ 
eral  (Ripley,  1981).  Lawson  required  that  the  smallest  of  the  six  angles  in  the  two  tri¬ 
angles  larger  for  this  division  of  a  convex  quadrilateral  than  that  given  by  the  other 
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diagonal,  as  shown  in  figure  80.  The  left-hand  triangulation  is  chosen  because  angle  A 
is  larger  than  angle  B.  Thus  when  the  four  data  measurement  locations  are  not  used  in 
the  model  development,  then  the  triangulation  is  considerably  changed  and  may  not  be 
optimum.  Further,  Akima  (1984)  has  shown  that  poor  estimates  of  partial  derivatives 
usually  occur  when  a  thin  (or  slim)  triangle  is  involved  in  the  interpolation,  which  af¬ 
fects  Ae  accuracy  of  the  interpolation.  This  is  more  likely  to  have  occurred  when  the 
four  data  locations,  particularly  the  two  high  latitude  sites,  are  left  out. 
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Figure  80.  The  Lawson  criteria.  The  left-hand  triangulation  is  chosen  because 
angle  A  is  greater  than  B. 

Several  of  the  contour  plots  giving  the  correction  factors  in  flgures  7  through  30 
have  values  larger  or  smaller  than  the  extreme  values  used  in  the  interpolation  to  ob¬ 
tain  these  figures.  The  flgures  with  the  minimum  extreme  less  than  the  input  minimum 
extreme  usually  had  this  extreme  occurring  at  one  or  two  of  the  same  locations.  These 
two  locations  are  either  in  the  very  most  Eastern  part  of  Asia  or  in  the  North  Atlantic 
just  south  of  Greenland.  These  extremes  likely  would  not  have  occurred  had  the  data 
at  Thule,  Greenland,  been  used;  in  this  case  the  triangulation  would  have  been  consid¬ 
erably  different.  The  flgures  with  the  maximum  extreme  greater  than  input  maximum 
extreme  usually  had  this  extreme  occurring  at  one  or  two  of  the  same  locations.  The 
flrst  location  was  usually  in  the  vicinity  of  Guam.  The  second  location  was  in  the 
South  Atlantic  between  South  America  and  Africa.  This  latter  location  would  have 
been  affected  by  the  triangulation  if  the  data  from  Ibadan,  Nigeria,  and  Byrd  Station, 
Antarctia,  had  been  included.  However,  the  extreme  occurring  in  the  vicinity  of  Guam 
is  due  to  the  gradients  occurring  in  the  input  data  in  eastern  Asia.  The  data  there  goes 
from  large  negative  corrections  in  the  north  to  large  positive  corrections  at  Singapore. 
The  interpolation  program  extrapolated  to  obtain  the  large  positive  correction  factors  in 
the  vicinity  of  Guam.  Further,  in  this  case  the  triangulation  would  not  have  changed  if 
the  data  for  the  other  four  locations  had  been  included  in  the  model  development. 
Whether  this  correction  is  valid  or  not  is  unclear.  Additional  input  data  in  the  vicinity 
of  Guam  would  have  been  useful. 
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5.0  PROPOSED  PROCESS  TO  DEVELOP  A  CORRECTED 
I  MHz  RADIO  NOISE  MODEL 

A  three-step  process  is  proposed  to  develop  a  new  1-MHz  noise  model.  The  first 
step  is  to  obtain  correction  factors  for  additional  locations  to  increase  the  accuracy  of 
the  interpolation.  This  includes  adding  data  measured  at  Thule,  Greenland;  Byrd  Sta¬ 
tion,  Antarctica;  Ibadan,.  Nigeria;  and  Bill,  Wyoming.  Attention  will  be  given  to  the 
elimination  of  man-made  noise  at  these  locations.  The  second  step  is  the  interpolation 
of  the  data  to  a  100-latitude  by  84-longitude  grid  for  each  time  block/season  for  nu¬ 
merical  mapping.  Consideration  will  be  given  to  using  the  interpolation  method  used  by 
Spaulding  and  Washburn  (1985).  However,  a  two-dimensional  method  of  interpolation 
(Akima,  1978a,  1978b,  1979)  and  more  recent  algorithms  (Cline  &  Renka,  1984;  Renka 
1984a,  1984b,  1984c;  Renka  &  Cline,  1984)  will  be  examined  for  possible  use.  The  fi¬ 
nal  step  is  the  numerical  mapping  of  the  data  itself  to  produce  the  final  model.  A  nu¬ 
merical  mapping  technique  from  Zacharisen  and  Jones  (1970)  will  be  used.  However, 
the  mapping  will  be  done  in  local  time  rather  than  in  universal  time  as  applied  by 
Zacharisen  and  Jones. 

5.1  ADDITIONAL  DATA  LOCATIONS 

There  are  nine  possible  locations  for  which  additional  data  might  be  used  to  in¬ 
crease  the  accuracy  of  the  1-MHz  atmospheric  noise  model.  The  first  four  are  the  four 
original  1-MHz  measurement  locations  not  used  to  obtain  correction  factors.  These  are 
Thule,  Greenland;  Byrd  Station,  Antarctica;  Ibadan,  Nigeria;  and  Bill,  Wyoming.  For 
Thule  and  Byrd  Station,  the  effects  of  man-made  noise  will  be  removed  by  techniques 
described  below.  For  Ibadan,  a  correction  factor  of  zero  will  be  used  since  there  is  no 
new  data  available  for  this  site  from  that  already  used  in  the  original  CCIR  322  model. 
For  Bill,  the  correction  factors  for  Boulder  will  be  used  as  Spaudling  and  Washburn 
indicate  that  they  differ  little  from  those  at  Bill.  If  data  from  the  four  Soviet  Union  lo¬ 
cations  (Simferopol,  Sverdlovsk,  Tbilisi,  and  Kiev)  not  used  in  the  development  of 
CCIR  Report  322-3  can  be  used,  the  accuracy  in  the  northern  latitudes  would  be  in¬ 
creased.  These  data  will  be  reexamined  for  that  purpose.  The  latter  three  sites  are  lo¬ 
cated  close  to  either  Moscow  or  Ashkabad  and  may  not  provide  as  much  improvement 
as  would  Sverdlovsk.  One  additional  site  is  at  Ping-Cheng,  Taiwan,  China,  about  33  km 
from  Taipei.  Atmospheric  noise  measurements  have  been  made  at  this  location  using 
the  ARN-2  atmospheric  noise  measurement  equipment  from  1967  through  1973 
(Huang,  1977).  If  that  data  can  be  obtained,  then  the  need  for  data  in  the  vicinity  of 
Guam  may  not  be  necessary.  The  locations  of  these  nine  sites  are  noted  in  table  10. 
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Table  10.  Additional  atmospheric  noise  measurement  locations. 


Bill,  Wyoming 

105.2W, 

43.2N 

Byrd  Station,  Antarctia 

120.0W, 

80.0S 

Ibadan,  Nigeria 

3.9E, 

7.4N 

Kiev,  USSR 

30.3E, 

50.72N 

Ping-Cheng,  Taiwan 

121. 23E, 

24.95N 

Simferopol,  USSR 

34.03E, 

45.02N 

Sverdlovsk,  USSR 

61.07E, 

56.73N 

Tbilisi,  USSR 

40.0E, 

41.72N 

Thule,  Greenland 

68.7W, 

76.6N 

Data  measured  on  board  the  ship  USNS  Eltanin  might  be  of  use  in  verifying  the 
accuracy  of  the  interpolation  method  chosen  to  interpolate  the  correction  factors  to  a 
100-longitude  by  84-latitude  grid.  This  ship  made  measurements  from  April  1962 
through  November  1964  at  various  locations  off  the  coast  of  South  America.  The  ship 
usually  made  measurements  at  5  locations  during  a  month,  repeating  some  locations 
during  another  month  of  the  same  season.  As  a  result,  more  than  a  dozen  locations 
were  used  to  make  measurements. 

Atmospheric  noise  data  being  collected  by  Stanford  University  might  be  of  use  in 
verifying  ^e  accuracy  of  the  atmospheric  noise  model  at  VLF  frequencies.  Stanford 
University  is  presently  operating  a  global  network  of  eight  computer-controlled  receiv¬ 
ing  systems  for  the  measurement  of  radio  noise  in  the  10  to  32,000  Hz  (ELFATP)  fre¬ 
quency  band  (Fraser-Smith  et  al.,  1987).  Data  are  being  recorded  on  16  frequencies, 
including  10.2  and  32  kHz.  The  locations  of  the  sites  are  given  in  table  11  and  figure 
81.  Several  of  these  locations  are  in  areas  for  which  no  other  data  exists. 

Table  11.  Geographical  coordinates  of  ELFAT-F  noise  measurement  sites. 


Station 

Coordinates 

#• 

Thule,  Greenland 

77»N,  69'‘W 

Sdndre  Strdmfjord,  Greenland 

67»N,  51»W 

« 

New  Hampshire 

44»N,  72»W 

L'Aquilla,  Italy 

42»N,  13»E 

Stanford,  California 

37»N,  122“W 

Kochi,  Japan 

33-^  133‘'E 

Dunedin,  New  Zealand 

46»S,  170»E 

Arrival  Heights,  Antarctica 

78»S,  167»E 
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Hgure  81.  Locations  of  the  eight  ELFA^  radiometer  locations  on  a  world  map. 


5.2  REMOVAL  OF  THE  EFFECTS  OF  MAN-MADE  NOISE 

Basically,  three  methods  are  available  for  processing  the  measured  data  to  remove 
the  effects  of  man-made  noise  in  the  determination  of  the  noise  at  1  MHz  at  each 
measurement  location.  These  methods  can  be  used  together  or  alone  depending  on  the 
particular  case.  Having  obtained  the  noise  value  at  1  MHz,  the  correction  factors  are 
determined  as  described  in  section  3.3. 

5.2.1  Extrapolation  up  to  1  MHz 

This  method  was  described  in  section  3.2.  The  measured  noise  values  for  each  time 
block  and  season  are  graphed  as  a  function  of  frequency.  Then  the  atmospheric  noise 
frequency-dependence  curves  for  this  same  time  block/season  and  the  noise  curves  for 
quiet  rural  and  rural  man-made  noise  are  superimposed  on  this  curve.  By  examining 
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the  results,  it  can  be  determined  whether  the  measured  values  at  the  HF  frequencies 
are  contaminated  by  man-made  noise.  This  can  be  determined  by  comparing  the  slope 
of  the  values  at  HF  to  that  of  the  man-made  curves  and  by  comparing  the  magnitude 
of  the  noise  relative  to  that  of  the  man-made  curves.  If  it  appears  that  there  is  no  con¬ 
tamination,  then  the  frequency-dependence  curves  can  be  used  to  interpolate  the  data 
to  1  MHz.  If  it  appears  that  there  is  contamination,  then  the  data  points  below  1  MHz 
can  be  used  to  extrapolate  up  to  1  MHz.  Only  those  data  points  not  having  a  slope  the 
same  as  the  man-made  curve  are  used.  The  frequency-dependence  curve  that  fits  the 
data  points  is  determined.  Then  this  curve  is  used  to  extrapolate  up  to  1  MHz.  Exam¬ 
ples  are  given  in  figures  5  and  6,  which  are  discussed  in  section  3.2. 

A  graphics  package  such  as  Axum  for  personal  computers  could  be  used  to  real 
advantage  for  this  purpose.  Each  of  the  sets  of  frequency-dependence  curves  for  both 
atmospheric  and  man-made  noise  for  each  time  block/season  can  be  stored  as  an  im¬ 
age  graph  to  be  used  later.  Then,  after  the  data  for  a  site  for  a  time  block/season  have 
been  graphed,  the  appropiate  image  graph  can  be  superimposed.  Having  determined 
the  particular  frequency-dependence  curve  to  use  to  extrapolate  to  1  MHz,  Axum’s  ca¬ 
pability  to  zoom  in  on  a  portion  of  a  graph  while  in  the  “VIEW”  mode  can  be  used  to 
determine  an  accurate  estimate  of  the  1-MHz  value. 

The  method  described  here  has  the  disadvantage  of  only  using  the  frequencies  be¬ 
low  1  MHz.  Thus,  in  most  cases  half  the  data  set  is  lost.  If  this  method  is  used  with 
the  methods  described  below,  it  may  be  possible  to  use  the  data  above  1  MHz. 

5.2.2  Use  of  and 

If  the  parameters  and  Lj  are  available  for  a  particular  measurement  value  of 
noise,  they  can  be  used  to  determine  if  the  measured  noise  value  is  contaminated  by 
man-made  noise.  If  the  procedure  described  below  is  used,  and  it  is  determined  that 
the  measurement  is  contaminated,  then  the  measured  value  should  be  eliminated  from 
the  set  of  data  unless  some  other  procedure  is  used  to  remove  the  effects  of  man-made 
noise. 

The  measured  noise  parameter  that  will  reflect  any  contamination  by  man-made 
noise  will  be  the  logarithmic  parameter,  Le  (Crichlow  et  al.,  1960b).  This  contamina¬ 
tion  generally  will  cause  the  value  of  Lj  to  be  less  than  it  would  have  been,  had  the 
recorded  value  been  only  atmospheric  noise.  In  determining  the  amplitude-probability 
distribution  (APD)  from  the  three  measured  statistical  moments,  F„  Vj,  and  Ld 
(Crichlow,  Disney,  &  Jenkins,  1960a;  Crichlow  et  al.,  1960c),  contaminated  values  of 
Ld  may  be  found  that  will  not  give  a  solution  of  the  APD.  That  is,  for  a  given 
only  a  certain  range  of  values  of  Ld  will  result  in  an  APD  representative  of  atmos¬ 
pheric  noise.  The  combination  of  a  given  Va  and  a  minimum  allowable  Ld  results  in  a 
nonunique  solution  (i.e.,  there  is  an  infinite  number  of  solutions  all  with  the  same  Va 
and  Ld  combination).  In  figure  82  the  lower  curve  on  the  graph  indicates  the  minimum 
value  of  Ld  that  will  give  an  APD  by  the  method  of  Crichlow  et  al.  (1960a)  for  a  given 
value  of  Vd.  It  can  therefore  be  used  to  determine  if  the  measured  noise  value  is 
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contaminated.  The  upper  curve  in  the  graph  is  a  curve  of  the  most  probable  value  of 
Ld  for  a  given  Vj  determined  as  the  best  fit  for  the  integrated  moments  from  over  60 
measured  APDs  of  uncontaminated  atmospheric  radio  noise. 

5.2.3  Estimates  of  Man-Made  and  Galactic  Noise 

If  estimates  or  measurements  of  man-made  noise  and  galactic  noise  are  available,  it 
should  be  possible  to  remove  the  effects  of  these  noise  sources  by  using  the  inverse  of 
the  normal  procedure  used  to  combine  the  estimates  of  noise  from  the  three  main 
sources  of  noise  used  in  prediction  programs.  In  the  case  of  Byrd  Station,  Enkoping, 
and  Thule  measurement  locations,  Herman  (1962,  1963,  1964)  has  determined  the 
man-made  noise  environment.  For  galactic  noise,  which  is  only  received  at  frequencies 
above  the  vertical  critical  frequency  of  the  F-region,  it  is  necessary  to  know  the  critical 
frequency  of  the  F-region  for  the  hour  and  month  of  the  measurement.  F-region  critical 
frequency  data  are  available  on  computer  media  for  this  purpose  for  the  following 
noise  measuremnt  sites:  Alma  Alta,  Ashkabad,  Boulder,  Byrd  Station,  Ibadan,  Irkutsk, 
Khabarovsk,  Moscow,  Murmansk,  New  Delhi,  Ping-Cheng,  Pretoria,  Rabat,  Simferopol, 
Singapore,  Sverdlovsk,  Tbilisi,  and  Thule.  This  includes  all  but  two  entries  in  table  10. 

5.2.4  Use  of  High-Latitude  Absorption  Events  to  Determine  Man-Made  High-Latitude 
Man-Made  Noise  Environments 

In  both  northern  and  southern  Arctic  regions,  Herman  (1962,  1963,  1964)  has 
shown  that  atmospheric  noise  decreases  during  polar  cap  absorption  (PCA)  events, 
most  markedly  in  the  HF  band.  He  has  estimated  the  probable  magnitude  of  man¬ 
made  noise  in  the  HF  band  from  data  taken  when  atmospheric  noise  is  absent.  He  es¬ 
timated  the  man-made  noise  level  during  PCA  events  at  both  Thule,  Greenland,  and  at 
Byrd  Station,  Antarctica. 
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Herman  used  a  noise  index,  defined  as  the  daily  average  fractional  departure  of  the 
measured  noise  from  the  month-hour  median,  to  select  anomalous  absorption  periods 
between  March  1958  and  December  1959.  Both  polar  cap  and  auroral  zone  absorption 
events  were  detected  with  this  index.  At  Byrd  Station,  14  PCA  events  were  identified 
that  correspond  to  previously  reported  Nothem  Hemisphere  PCAs.  It  was  shown  that 
PCA  has  a  remarkable  effect  on  HF  atmospheric  radio  noise.  The  noise  level  is  se¬ 
verely  depressed, .especially  on  5.0  and  10  MHz,  and  remains  so  during  the  entire  PCA 
event.  The  noise  level  is  depressed  as  much  as  25  dB  below  the  median  during  PCA. 


Figure  83  shows  the  mimimum  noise  on  each  frequency  during  all  PCAs  for  Byrd 
Station  averaged  together  and  plotted  as  a  function  of  frequency.  (The  minimum  noise 
recorded  during  absorption  events  in  all  cases  was  the  the  minimum  recorded  in  the 
month  without  regard  to  ionospheric  conditions.)  The  average  deviation  is  ±3  dB  or 
less  on  all  frequencies.  Curve  B  shows  the  minimum  radio  noise  level  that  can  be  ex¬ 
pected  at  any  time  at  Byrd  Station,  regardless  of  the  source.  Because  ionospheric  ab¬ 
sorption  is  maximum  at  the  gyro-frequency,  Herman  felt  that  the  measured  noise 
power  on  2.5  and  5.0  MHz  during  intense  absorption  events  consisted  entirely  of  man¬ 
made  noise.  Therefore,  he  estimated  a  man-made  noise  curve  (curve  C),  with  the  same 
slope  as  used  in  CCIR  for  man-made  noise.  At  10  and  20  MHz,  the  noise  power  is 
above  the  indicated  man-made  noise  level  because  the  absorption  of  the  distant  atmos¬ 
pheric  noise  changes  inversely  as  the  square  of  the  frequency  and  because  of  a  large 
contribution  from  cosmic  noise  at  frequencies  near  and  above  the  F-layer  critical  fre¬ 
quency.  Below  2.5  MHz,  the  noise  power  becomes  increasingly  greater  than  the  man¬ 
made  level  because  complete  absorption  of  the  atmospheric  noise  does  not  take  place 
on  the  frequencies  that  do  not  penetrate  completely  throughout  the  D-region  before  re¬ 
flection.  Curve  A  is  the  median  autumn  00-04  time  block  curve  for  Byrd  Station  for 
comparison. 


Of  all  the  atmospheric  noise  measurement  sites,  only  Murmansk  is  far  enough 
north  to  experience  a  PCA  event.  The  geomagnetic  latitude  there  is  64.1  degrees  north. 
This  means  reprocessing  the  noise  data  for  this  location,  but  likely  will  improve  the 
values  of  1-MHz  atmospheric  noise  estimates  for  this  site.  It  does  require  the  daily 
noise  measurements  to  determine  the  noise  index  to  establish  the  days  for  which  PCA 
events  occured  while  noise  measurements  were  underway. 


5.3  INTERPOLATION  OF  CORRECTION  FACTORS  TO  THE  STANDARD  GRID 


After  obtaining  the  correction  factors,  the  next  step  is  to  interpolate  these  data  to  a 
standard  100-latitude  by  84-longitude  grid.  Once  this  has  been  accomplished,  then  the 
correction  factors  are  added  to  noise  values  for  the  original  CCIR  Report  322  also 
given  at  these  same  grid  locations  for  each  time  block  and  season. 
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Figure  83.  Minitnum  noise  level  observed  at  Byrd  Station  during  PCA. 


To  interpolate  the  data  to  a  standard  100-ladtude  by  84-longitude  grid  it  is  neces¬ 
sary  to  choose  an  accurate  interpolation  method.  The  available  methods  seem  to  fall 
into  two  categories.  The  first  category  includes  methods  that  interpolate  data  in  a 
plane.  Methods  that  could  be  used  that  fall  into  this  category  include  Akima  (1978a, 
1978b,  1979),  Cline  and  Renka  (1984),  Renka  (1984c),  and  Renka  and  Cline  (1984). 
Additional  methods  are  described  by  Franke  (1982).  The  disadvantage  to  these  meth¬ 
ods  is  that  the  interpolation  is  affected  by  the  boundary  imposed  by  the  available  data 
points.  This  affects  both  the  number  of  triangles  and  the  choices  themselves  of  the  tri¬ 
angles  chosen  for  the  interpolation.  The  second  category  eliminates  the  boundary  by 
interpolating  the  data  on  a  sphere.  There  are  two  methods  available  for  this  purpose 
(Lawson,  1982;  Lawson,  1984;  Renka,  1984a;  Renka,  1984b).  The  method  due  to 
Renka  is  available  from  the  Association  for  Computing  Machinery.  As  the  atmospheric 
noise  data  was  measured  on  a  sphere,  these  two  methods  would  be  applicable.  It  will 
be  necessary  to  test  the  accuracy  of  these  two  methods  against  a  benchmark.  One  such 
benchmark  might  be  computed  results  from  a  model  of  the  Earth’s  magnetic  field. 
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5.4  NUMERICAL  MAPPING  OF  THE  NOISE  DATA 


The  final  step  in  the  development  of  an  improved  atmospheric  noise  model  is  the 
numerical  mapping  of  the  lOO-by-84  grid  data  sets  for  use  in  computer  programs.  Two 
numerical  mapping  techniques  are  available  for  this  purpose  (Lucas  and  Harper,  1965; 
Zacharisen  and  Jones,  1970).  The  Zacharisen  and  Jones  method  results  in  a  consider¬ 
able  reduction  of  coefficients  but  was  originally  designed  to  be  used  with  data  given  as 
a  function  of  universal  time.  The  reduction  in  coefficients  significantly  reduces  the 
computation  time  of  models  produced  by  the  technique.  Sailors  (1985)  found  that  the 
local-time  Lucas  and  Harper  model  was  not  always  more  accurate  than  the  universal¬ 
time  Zacharisen  and  Jones  model.  Moreover,  Sailors  reported  the  development  of  a 
1-MHz  minicomputer  atmospheric  noise  model  in  local  time  using  a  modification  of 
the  Zacharisen  method.  He  found  that  for  each  and  every  case,  the  local-time  map  had 
a  lower  error  than  the  corresponding  universal-time  map.  This  occurs  because  the  gra¬ 
dients  of  atmospheric  noise  as  a  function  of  longitude  are  much  smaller  than  those  in 
universal  time.  Hence,  it  appears  that  the  Zacharisen  and  Jones  numerical  mapping 
technique  applied  in  local  time  is  more  accurate  than  the  Lucas  and  Harper  method 
and  would  be  the  method  to  use. 

One  major  problem  exists  with  numerical  mapping  techniques.  That  is  the  tendency 
to  smooth  out  physical  properties  of  the  atmospheric  noise,  particularly  at  low  latitude 
where  the  gradient  in  atmospheric  noise  maximizes  due  to  concentration  of  noise 
sources  in  that  region.  As  the  number  of  harmonics  in  the  numerical  map  is  reduced 
to  make  the  model  efficient,  the  high-level  contours  begin  to  disappear,  enhancing  the 
problem.  Sailors  (1983)  showed  that  an  improvement  can  be  made  by  making  a  trans¬ 
formation  so  that  the  latitude  input  into  the  latitude  function  in  the  numerical  mapping 
process  is  spread  further  away  from  the  equator  than  the  actual  latitude.  This  transfor¬ 
mation  resulted  in  a  more  uniform  variation  of  atmospheric  noise  as  a  function  of  the 
modified  latitude  than  for  the  actual  latitude.  The  transformations  used  preserved  the 
stability  of  the  latitude  function  at  the  poles.  One  problem  was  that  the  transformations 
were  applied  with  the  universal-time  version  of  the  Zacharisen  and  Jones  mapping 
technique.  The  gradients  in  longitude  had  a  tuidency  to  introduce  errors  at  other  longi¬ 
tudes  removed  from  the  atmospheric  noise  storm  centers.  When  used  with  the  local- 
wime  Zacharisen  and  Jones  technique,  it  can  be  expected  to  work  even  better  in  im¬ 
proving  the  atmospheric  noise  model. 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  has  presented  the  probable  cause  and  a  recommended  solution  for  a 
discrepancy  in  the  CCIR  Report  322-3  radio  noise  model.  This  report  has  not  deter¬ 
mined  whether  there  is  a  discrepancy  between  the  new  CCIR  model  and  the  data  val¬ 
ues  used  to  develop  it.  Nor  has  there  been  an  attempt  to  determine  the  validity  of  the 
measured  data  values  used  to  develop  the  model  or  to  validate  the  model  against  any 
other  data. 

The  basis  for  this  discrepancy  was  found  to  be  in  the  procedure  used  to  prepare 
the  measured  noise  data  for  the  determination  of  a  global  numerical  representation  of 
the  1-MHz  data.  The  procedure  followed  in  the  development  of  the  model  was  to  de¬ 
termine  correction  factors  to  the  old  CCIR  model  for  each  measurement  site,  to  inter¬ 
polate  these  corrections  to  a  100-latitude  by  84-longitude  grid  for  each  time  block/sea¬ 
son,  to  add  the  correction  factors  at  each  grid  point  to  corresponding  values  for  the  old 
CCIR  model,  and  finally  to  numerically  map  the  resulting  data  for  each  time  block  and 
season.  Nineteen  locations  were  used  in  the  final  model.  Four  sites  used  in  the  original 
CCIR  model  were  not  used.  These  include  Bill,  Wyoming;  Byrd  Station,  Antarctica; 
Ibadan,  Nigeria;  and  Thule,  Greenland.  As  no  correction  factors  were  obtained  for 
these  locations  or  a  correction  factor  of  zero  was  used,  the  interpolation  algorithm  used 
to  obtain  the  100-latitude  by  84-longitude  grid  of  correction  factors  supplied  other  val¬ 
ues.  For  Bill,  Wyoming,  the  result  is  not  too  serious,  but  for  the  other  three  sites  the 
error  is  at  some  seasons  and  time  of  day  serious.  For  Thule,  Greenland,  the  maximum 
and  minimum  errors  in  the  correction  contours  were  10.1  and  -10.8  dB,  respectively. 
For  Ibadan,  Nigeria,  the  maximum  and  minimum  errors  were  12.5  and  -1.5  dB,  re¬ 
spectively,  For  Byrd  Station,  Antarctica,  the  maximum  and  minimum  errors  were  12.0 
and  3.0  dB,  respectively.  Examination  of  the  geographical  extent  of  these  errors  re¬ 
veals  that  the  error  is  not  confined  to  the  measurement  location  but  in  fact  is  very 
large.  It  was  found  that  the  error  as  a  function  of  frequency  was  diumally  dependent. 
An  error  of  10  dB  at  1  MHz  was  more  serious  at  another  frequency  during  local  day¬ 
time  than  at  night.  Finally,  the  absence  of  the  data  locations  affected  the  accuracy  of 
the  interpolation  itself. 

Because  of  the  errors  in  the  CCIR  Report  322-3  atmospheric  noise  model,  it  is  rec¬ 
ommended  that  it  be  used  with  caution.  It  is  most  accurate  in  Europe,  Asia,  the  Indian 
Ocean,  the  western  Pacific  from  Asia  to  the  date  line,  and  Australia.  It  is  most  inaccu¬ 
rate  in  both  the  northern  and  southern  high  latitudes,  the  Arabian  Peninsula,  northern 
Africa,  and  the  mid-Atlantic  Ocean  area.  For  applications  in  the  latter  list  of  areas,  the 
user  should  consider  using  the  original  CQR  Report  322  model. 

A  three-step  process  was  proposed  to  develop  a  new  1-MHz  atmospheric  noise 
model.  The  first  step  was  to  obtain  correction  factors  for  additional  locations  to  in¬ 
crease  the  accuracy  of  the  interpolation.  Nine  additional  locations  were  identified  as 
locations  for  which  there  might  be  sufficient  data  to  obtain  correction  factors. 

Three  methods  were  presented  that  might  be  used  in  removing  the  effects  of  man¬ 
made  noise  from  the  noise  data.  Second,  two  alternate  methods  of  interpolation  were 
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identified.  It  was  proposed  that  the  accuracy  of  these  methods  be  examined  against  a 
benchmark  to  be  determined.  Finally,  it  was  proposed  that  the  local  time  version  of  the 
Zacharisen  and  Jones  (1970)  numerical  mapping  technique  be  used  to  develop  the  final 
model.  In  addition,  it  was  proposed  that  a  latitude  transformation  be  used  to  increase 
the  accuracy  of  the  numerical  mapping  itself. 

It  is  recommended  that  the  U.S.  counterpart  to  CCIR  Study  Group  6  take  a  lead 
role  in  submitting  a  corrected  model  to  the  CCIR. 
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ment  location  but  in  fact  is  veiy  large.  It  was  found  that  the  error  as  a  function  of  frequency  was  diurnally  dependent.  An  error  of 
10  dB  at  1  MHz  was  more  serious  at  another  frequency  during  local  dqytime  than  at  night.  The  absence  of  the  data  locations  affected 
the  accuracy  of  the  interpolation  itself. 

The  recommended  course  of  action  includes  (1)  obtain  correction  factors  for  additional  locations  to  increase  the  accuracy  of  the 
interpolation;  (2)  test  the  method  of  interpolation  against  a  benchmark;  (3)  use  the  Zacharisen  and  Jones  numerical  mapping  tech¬ 
nique  applied  in  local  time  to  develop  a  new  1-MHz  model;  (4)  consider  using  a  latitude  transformation  to  increase  the  accuracy  of 
the  numerical  mapping  technique;  (5)  submit  a  corrected  model  to  the  CCIR;  and  (6)  use  the  CCIR  Report 322-3  atmospheric  model 
with  caution,  especially  in  the  northern  and  southern  high  latitudes,  the  Arabian  Peninsula,  northern  Africa,  and  the  mid-Atlantic 
Ocean  area.  In  these  cases,  consider  using  CCIR  Report  322  instead  until  an  improved  model  is  available. 
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